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Abstract 

In our paper we show that the Camassa-Holm equation does not represent a long wave 
asymptotic due to a major inconsistency with the theory of shallow water waves. We state that 
any solution of the Camassa-Holm equation, which is not asymptotically close to a solution 
of the Korteweg-de Vries equation is an artefact of the model and irrelevant to the theory of 
shallow water waves. 
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1 Introduction 

A partial differential equation 



2toU y + U T + 3UyU - 
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known in the literature as Camassa-Holm equation is a fine example of an integrable system with 
many interesting and rather unusual properties. As an i ntegrable equation it was discovered in 
Fokas fe Fuchssteiner ( 1980h . Fokas fe Fuchssteinerl dl981 ). but equation {J) had not been really 
noticed until the publication ICamassa fc Holml ( 19931 ) where the authors proposed it as a model 
for shallow water waves. Moreover, Camassa and Holm have shown that in the case u> — equation 
(HJ) possesses a peculiar solution with a cusp 



U = i>exp(— \y — vt\) 



(2) 



which they called a "peakon" . Moreover, they also found exact multi-peakon solutions. The paper 
of Camassa and Holm has triggered an enormous avalanche of publications with the developments 
of mathematical theory for this new class of integrable systems and speculations about possible 
applications of the Camassa-Holm equa tion to the theory of shallow water waves (including the 
problem of tsunami. ILakshmanan ( 20071 )) and other long wave asymptotic theories. 

There were a number of publications where authors raised some criticism to the original deriva- 
tion and proposed new versions of the derivations starting from the reduction of the Green-Naghdi 
model, reduction of the generalised Serre model or a direct multis c aling a symptotic expansion of 
the E ul er equation (see f or example Johnson ( 2002 V_ Dullin et al. ( 2003 ). Constantine fc Lanned 
(2009), iDias fc Milewskil (|2010l) ICamassa et all (119941) 1. The main concern of these derivations 
was to achieve a fine adjustment of the coefficients in the first few terms of the asymptotic ex- 
pansion in order to make the equation corresponding to the truncated expansion integrable. The 
aim of this paper is to point out that the Camassa-Holm equation does not correspond to any 
dominant balance in the asymptotic expansion and thus it does not represent a long wave asymp- 
totic for water waves. We have analysed the above derivations and found them inconsistent with 
the basic principles of asymptotic theory. Indeed, one cannot keep the first principal asymptotic 
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contribution (corresponding to the Korteweg-de Vries theory) together with the next asymptotic 
correction, then truncate the expansion and balance these two contributions by a re-scaling. Such 
a re-scaling contains the small asymptotic parameter and thus violates the original assumption 
about the ratio of the water depth and the characteristic wave length. We have computed the next 
(neglected in the above derivations) term of the expansion and have shown that after the proposed 
re-scaling it is of the same order as any term accounted in the Camassa-Holm equation. Moreover, 
the parameter u in ([T]) appears in denominators of the neglected terms and thus cannot be set to 
zero (for the existence of peakon solutions) in any long-wave asymptotic theory. Here we should 
mention that the fact that the pea kon solut i ons ar e irrelevant fo r shall ow water waves was well 
understood earlier and published in IJohnsonl (|2002t ). iDullin et all (|2003l ). Using the exact soliton 
solution of the Camassa-Holm equation we have shown that the neglected terms are of the same 
order as terms accounted in equation (JT]). 



2 Long wave expansion and the Camassa-Holm Equation 



In this section, following IWhithaml t|l974h . we give a sketch of long wave asymptotic expansion be- 
yond the Korteweg -de Vries (KdV) th e ory . We shall illustr ate the derivation of the Camassa-Holm 
equation following IDullin et all (|2003l ). IDullin et all (2004) and point out where the inconsistency 
occurs. Also we will show that the terms neglected in the theory of the Camassa-Holm equation 
are of the same order as any term accounted in ([1]). We claim that the inconsistency cannot be 
removed if one uses a reduction of the Green-Naghdi model (whi ch itself i s an a pproximation), or 
by choosing the value of the velocity potential inside of the flow ()johnson ( 20021 )1 or by any other 
method. We shall neglect the surface tension, since it does not affect our argument but simply 
makes expressions look more complicated. Also we shall neglect viscosity and compressibility of 
water. 

Two dimensional motion of vorticity free fluid is described by the potential <f>(x' ,z',t') of the 
velocity field u(a:', z', t') = V$. In the bulk of the fluid the potential satisfies the Laplace equation 



®x>x> +<£z>z> =0 

with the boundary condition $ z / — at the bottom z' — 0. At the free surface of the fluid 
z' = ho + H(x',t') there are kinematic and dynamic boundary conditions 



where g is the acceleration due to gravity and Jiq is the undisturbed depth of water. 

Long wave asymptotic expansion assumes a small parameter e = 1iq/L 2 where L is a typical 
wavelength of the wave. Another dimcnsionlcss small parameter of the theory is /i = ao/tiQ where 
a is a typical amplitude of the wave. We shall assume that e<n« l . Actually one can s e t /x = e 
and develop the theor y with one parameter, but following iJohnson (|2002l) . IDullin et all (120031) . 
Constantine fc Lannes! (2009) we shall keep both parameters for better control over the terms. 

Introducing dimensionless variables 



x! = L.t, z' = llnZ, t' = 



-.t, H = fih T], $ = fil^/ghoc, 
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we re-write the above system of equations in the form: 



= 



= 



1 



Vi 
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(3) 
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Starting from here we shall develop asymptotic expansion as a series in e n fj, m , n,m G Z + . 
We shall illustrate the derivation the Camassa-Holm equation with corrections in three steps. We 
begin wi th the derivation of the Boussinesq expansion up to the order c n fj, m , n + m = 3. Then, 
following IWhithaml (|l974l ). we make a reduction to the KdV theory describing unidirectional wave 
propagation. Finally we transform the equation obtained to the form ([T]), keeping terms of order 
e n /z m , n + m = 3 to demonstrate that the Camassa-Holm equation does not represent a long wave 
asymptotic for surface waves. 



2.1 The Boussinesq expansion 

The Boussinesq expansion aims to eliminate the dependence on the vertical coordinate z and reduce 
Q - P to a system of equations on the elevation r\ = rj(x,t) and the horizontal component of 
the velocity field at the bottom w = 4>(x,0,t) x . In this Section we shall follow the construction 
presented in detail in IWhithaml (jl974 (Chapter 13.11), but will keep more terms in the expansion. 
It follows from the Laplace equation (|3|) and the boundary condition at the bottom ([4]) that 



(x,z,t) 
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(7) 



(2n)\ dx 2r 

where F(x, t) is the value of the potential 4> at the bottorrO], and thus w = F x . 

In order to reduce the system (|3 j) - ([B]) to two equations for functions r/(x,t) and w(x,t) we 
substitute cj)(x,z,t) (JTJ), in ([S]) and ©, then we differentiate in x the equations obtained from ^ 
and replace F x by w in the both equations. Keeping terms of order e"/i m , n + m < 3 we get 

= 7) t + w x + n (r)w) x ~ ^w xxx 
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(8) 
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4 w x w xxx + w xxxx w) - — - w t 
x 720 



et al. (2003) the authors use a different geometry, namely the bottom is set at z = — 1. Their solution 
does not satisfy the boundary condition at 2 = — 1 (presum ably due to a mispr int). To rectify the misprint one has 
to replace z by z + 1 in the right hand side of (2.7), (2.8) in lDullin eTai] j2003f) . 
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respectively. First two lines in ([H]),© coincide with equation (2.9) in lDullin et all (j200 3rl. For the 
purpose of our paper we are keeping the next order in the expansion. There are no obstructions 
to find higher order terms if required. 



2.2 Reduction To Unidirectional Waves. The KdV theory with higher 
asymptotic corrections 

System describes waves propagating in both direction. Th ere are many ways to reduce it 

to unidirectional wave p ropagation . In o rder to be consistent with iDullin et all (|2003l ) we employ 
the method proposed in lWhitham ( 1974 ). Namely, we shall assume 



fc=l n=0 



and request that equations ([I]),® coincide upon this assumption, that would enable us to determine 
the coefficients /&„ [rj] . Keeping terms with fc < 3 we get 
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Here D^ 1 denotes integration. Assuming rfc — > rapidly enough as x — > -co one c an se t 
Dx 1 (Vx) ~ f! o r lxdx. T he first line of the e xpan sion (HU1) one c a n find in IWhithaml (|1974l ). 
terms jin were derived in Marchant fe Smyth ( 1990l) and Johnson ( 2002 ). here we extend the 
expansion to the terms f^ n of order /_t n e m , n + m — 3. 

Substitution of (flU)) in either © or © leads to equation 
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The first line of this expansion is the st andard Korteweg-de Vries e quation, c orrect ions in the 
second line is the well known result (see Marchant fc Smyth ( 1990t ). Johnson ( 20031 )). For the 
purpose of our paper we proceed to the terms of order fi n e m , n + m = 3. 



2 In lDullin et all i2003T ) in the first equation (2.9) there is a misprint in the sign at the term proportional to <5 4 
(in our paper it is the term proportional to e 2 in J8t). 
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2.3 Asymptotic near-identity transformation 

Following IDuUin et al. (2003) we shall apply the Galilean and asymptotically invertible near- 
identity transformations to equation The purpose of these transformations is to bring the 
first two lines of equation (jlip in the form, which can be re-scaled to equation (JlJ. 
First we apply the Galilean transformation 

X = x-St, T = t, (12) 

where S is a constant which will be determined later. Then we perform the Kodama transformation 

r)(X,T) = u + fi {aiu 2 + a 2 u x D x 1 (u)) + e/3u X x (13) 

to a new dependent variable u = u(X,T). And finally we apply the Hclmholtz operator % = 
1 — ejd x to the equation obtained. 

In the Galilean transformation and the Helmholtz operator we set 6 — 9/19 and 7 = 19/60 in 
order to vanish the coefficients at the terms uxxx and uxxxxx respectively. The choice ai = 
7/20, 02 = —1/5, (3 = 1/30 excludes the term u 2 ux and guarantees (see details in Dullin et all 
(2003)) that the first line of the resulting equation 

10 3 fi 19 e fie . 

= u T + — u x + — uu x - utxx - 777: (38 u X xu x + 19 u X xxu) 
19 2 bO 120 

+ ^1200 Uxxxxxxx ~ UxxD x 1 («) (" 2 - 2 uxDx 1 (u)) 
+ t^tjtj (976 uuxuxx - ASuuxxxxDx 1 (u) + 48u XX 2 D x 1 (u) + 680u X xxu 2 

fie 2 

+2765 ux 3 ) + — — - (903 u X xxxu x + 316u X xxxxu + 305u X xu X xx) (14) 
' 3600 

can be re-scaled to (JTJ). 

Until this stage the asymptotic theory is consistent and e quation (1141) i s asym ptotically equiva- 
lent to the KdV expansion. The first line of (1141) was derived in Dullin et oil 42003(1 . our contribution 
is in the retaining of the next corrections in the asymptotic expansion (the last three lines in (j!4|l ). 



3 Derivation of the Camassa-Holm equation and its incon- 
sistency with the asymptotic expansion 

Analysing the publications with derivations of the Camassa-Holm equation as a long wave asymp- 
totic expansion we notice that they have a similar pattern. Starting from the Euler equation or 
a certain well established model of water waves (the Green-Naghdi model, the generalised Serre 
model, etc) the authors arrive to the equation similar (up to an inessential re-scaling with constant 
coefficients) to the first line of ([T4"|) . Then they truncate the expansion at this level and re-scale it 
to the form ([I]). We shall do the same re-scaling, but accounting the next correction. 

The most general re-scaling of variables that transforms the first line of equation (|14l) into (JTJ) 





AU, y=—^^X, t — —— _ T, (15) 



where A is an arbitrary constant and to — 10/(19 A fi). The re-scaling ([15]) balances the terms in 
the the equation by eliminating the small parameters e and fi (except the the term with oj). It is 
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easy to see that after this re-scaling the small parameters disappear from the correction (the last 
three lines in (fT4")l ). which takes the form 

12 9 fi 2440 

5^ U„U, (D- (U )f-—^U„D-^U) + — U,UU„ 

+ ^ {U - )3 + wz W on - ^ UU '- D '' «n + ^ &u„ 

1806 610 632 | 446 

y wyy ' ^w^vvv ' gg~[ ^^vvvvv > 2527 w yyyyyyy ' ' ' 

One can also demonstrate that the small parameters disappear from all higher corrections. 

One could hope that for the exact soliton or peakon solutions the correction term (fTl)|) and all 
higher corrections vanish. It is obviously not the case for peakons, moreover the constant u>, which 
has to be set zero (for the existence of peakons) is in the denominator. It does not hap pen with 
soliton s either. To show that we used the exact soliton solution of equation |T]) taken from Uohnsonl 



(120031) 



u (y> T ) = m: - ■ ; 2, 



1 + 2cjc- 1 sinh 2 (6») 

where c is an arbitrary constant satisfying condition c > 2ui and 9 is a function of y — ct implicitly 
given by equation 



29 ( cosh (9 — arctanh\/l — 2ljc 1 

y — ct = + In 1 



VT — 2ojc 1 \ cosh (0 + arctanhVl — 2cjc 1 ) I 

Taking c = 1, w = 0.2we compared a contribution from one of the terms of the Camassa-Holm 
equation with the value of the correction (fTB| . It is shown on Fig.l that the correction (|16l) is much 
bigger then a contribution from the term of the equation ([1]) evaluated on the soliton solution with 
this choice of parameters and this fact does not depend on the choice of the term. 

The re-scaling (p~5]) is a basic error, which leads to the inconsistency with the long wave asymp- 
totic theory. Indeed, suppose we have started from equation ([T]) and have found its solution of a 
size or characteristic wavelength A ~ 1. Re-scaling this solution to the variable X we find from 
(11511 that the size in this variable is A = -^§=\ 0.56^A. Since the Galilean transformation 

* 9 2V2S5 



(|T2j) docs not change the scale, the wave has the same size in the variable x. Coming back to the 
physical dimensional variable x' we realise that the size of the wave is A' = LA « 0.56A/io which is 
in contradiction with the long wave assumption /iq/(A') 2 ps 3.2/A 2 -C 1. 
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One can consider solutions of ([TJ of an extremely large characteristic length A 3> 1, but with 
the same accuracy such solutions can be described by the Korteweg-de Vries e quation with th e 
first correction ([Til (whos e integration is not much dif ferent from the KdV itself Kodamal (|l985 ) , 



Hiraoka fc Kodamal (|2009t )). For A > 1 it follows from lDullin et all (|2004l) that equations ([J) and 



(jlip are asymptotically equivalent, but the theory of the latter is much simpler and well developed. 



4 Conclusion 

The Camassa-Holm model for shallow water waves docs not represent the long-wave asymptotic. 
Any solution of this model, which is not asymptotically close to a solution of the Korteweg-de 
Vries equation is an artefact of the model and irrelevant to the theory of shallow water waves. In 
the literature there are many papers with implicit criticism of various aspects of the derivation 
of the Camassa-Holm model and its validity as well as attempts to rectify them, but they only 
contribute further to the confusion. Serious concerns abou t the asymptotic sense of the Camassa- 
Holm equation as a water wave theory has been raised in Johnson ( 2002h . In our paper we put 



an end to desperate attempts to justify the Camassa-Holm model using the long wave asymptotic 
theory. Based on our consideration it is not diffi cult to conclude, that neit her the Camassa-Holm 



model nor the Degasperis-Procesi equation (jDegasperis fc Procesil ill!)!)!); ; can represent a long 



wave asymptotic in problems of Hydrodynamic, Physics of Condensed Matter, Plasmas, etc, and 
that the peakon solutions are irrelevant for the long wave asymptotic theory. 

Having said so, we do not want to undermine the mathematical value of equation ([1]). It is a 
fine example of an integrable multi-Hamiltonian system with interesting associ ated spec t ral th eory, 
unu sual properties, etc. It ma y have other applications due to the discovery bv lMisiolekl (|l998h (see 



also Khesin fc Misiolek ( 20031) ) that equation JT]) is the Euler equation for the geodesic flow on the 



Virasoro group with respect to the right-invariant Sobolev H 1 -metric. As a mathematical object, 
the class of integrable equations discovered by Fuchssteiner and Fokas is a valuable contribution 
to the theory of differential equations. 
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